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Abstract 

The effect of molecular vibrations on electronic transport is investigated with the smallest 
fullerene C20 bridge, utilizing the Keldysh nonequilibrium Green's function techniques combined 
with the tight-binding molecular-dynamics method. Large discontinuous steps appear in the dif- 
ferential conductance when the applied bias-voltage matches particular vibrational energies. The 
magnitude of the step is found to vary considerably with the vibrational mode and to depend on 
the local electronic states besides the strength of electron-vibration coupling. On the basis of this 
finding, a novel way to control the molecular motion by adjusting the gate voltage is proposed. 

PACS numbers: 61.48. +c, 73.63.-b 
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Inelastic transport associated with local heating in nanoscale devices has been a growing 
interest in the fields of nanoscience and nanotechnology over the past few years 

mm 

1^. Fullerenes are considered promising candidates for basic elements in nanoscale devices, 
therefore, electronic transport in fuUerene bridges has received significant attention both 
e.pen.enta.ly flflQ and theoretically QQ. McE„e. and c„-wo..e. ,nea.n.ed cn„.ent- 

voltage (I-V) characteristics of the Ceo and C140 bridges 0, 3] and reported that coupling 
between electronic and vibrational degrees of freedom plays an important role in electronic 
transport. 

The strength of the electron-vibration coupling is known to be enhanced with decreasing 



the diameter of the fuUerene [jjj. Thus, among fullerenes, the smallest fuUerene C20 is 
expected to have the largest electron-vibration coupling. Recently, the C20 fuUerene was 
synthesized by Prinzbach et al To the best of our knowledge, the issue of inelastic 



transport through the C20 bridge in the presence of electron-vibration scattering has not 
been addressed. Although elastic transport through the C20 bridge has been theoretically 
investigated by several authors , characteristic features remain to be elucidated. 

The aim of this Letter is to clarify the role of molecular vibrations in electronic trans- 
port through the C20 connected to Au electrodes. A semi-infinite one-dimensional (ID) 
chain is used as a simple and ideal model for the Au electrode, as shown in Fig. ^ In 
an effort to achieve the aim, molecular vibrations in the C20 are first examined using the 
tight-binding molecular-dynamics (TBMD) method Q]. The influence of molecular vibra- 
tions on electronic transport characteristics of the C20 bridge is discussed in terms of the 
Keldysh nonequilibrium Green's function (NEGF) method Hf, a powerful method used to 
analyze inelastic transport in the presence of electron-vibration interactions 

Re- 
cently, the NEGF formalism has been successfully applied to inelastic transports in several 
nanostructures, such as atomic wires js^Ql and molecular junctions 13, Q, Hi- 

Before studying the vibrational properties of the C20 bridge, the stable structure of C20 
is determined using the TBMD method. The structure of C20 is distorted from the highest 
possible Ih symmetry to the D^d symmetry after the optimization, due to the Jahn- Teller 
effect. The optimized structure of the C20 contains three types of nonequivalent atoms, 
labeled a, 6, and c in Fig.^ There are four distinct bond lengths: aab =1.464 A, a^c =1.469 A, 
ttcc' =1.519 A, and acc" =1.435 A. The obtained bond lengths differ by less than 1 % from 
those obtained using the flrst-principles calculations 3|. 
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TABLE I: Vibrational energies (meV) of the isolated and connected C2o- 



Isolated C20 Connected C20 



Alg 


A2g 


Aiu 


A2u 




Eu 


A' 

^Ig 


A' 


58.2 


66.8 


69.2 


69.5 


29.9 


68.9 


62.6 


11.1 


98.6 


88.9 


74.0 


88.6 


56.1 


72.1 


100.8 


72.0 


105.2 


140.2 


131.5 


110.9 


76.8 


73.4 


106.6 


91.4 


135.1 




151.3 


127.6 


94.6 


92.4 


135.3 


112.4 


140.9 






148.8 


119.9 


113.3 


141.4 


127.7 


161.6 








132.6 


132.4 


161.6 


149.2 










137.7 


135.9 














148.8 


146.6 














163.6 


155.7 







The vibrational energies of the isolated C20 are first calculated by diagonalizing a 60x60 
force-constant matrix derived from the TBMD calculations. The results of the vibrational 
energies are classified using the irreducible representations of the D^cl symmetry group, as 
listed in Table HI Similarly, the vibrational energies of the C20 connected by two springs to 
the Au electrodes are calculated. The springs are attached to the a-atoms on either side of 
the principal axis of D^cL symmetry of the C20, as shown in Fig. The spring constant was 
assumed to be 4.37 eV/A^, an experimental value for the Ceo attached to Au electrodes Qj. 
In the present calculations, atomic vibrations of Au atoms were neglected because an Au 
atom is considerably heavier than a C atom. The connection of the electrodes only infiuences 
the vibrational modes with Aig and symmetries, because modes with other symmetries 
show no stretching of the springs between the C20 and Au electrodes. The vibrational 
energies shifted by the connection are listed in the A'^g and A'2^ columns in Table HI It is 
important to note that the vibrational mode of 11.1 meV in the A^^ column corresponds to 
the shuttle motion of C20 that goes back and forth between the two Au electrodes. Using 
these 55 vibrational modes, the Hamiltonian for the molecular vibrations of the C20 bridge 
can be written as =^ib = X]a^a(^a^a + where h\ {h\) is the creation (annihilation) 
operator of the vibrational quanta [i.e., phonon) with energy fvjj\. 
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The Hamiltonian for conduction electrons in the C20 bridge by the tight-binding model 
within the Hiickel approximation is subsequently described. The tight-binding Hamiltonian 
is expressed by the sum of five parts: = + '^ua + + '^mk + -^k- '^yi represents 
the Hamiltonian for an extended molecule, AU-C20-AU, including the edge Au atoms in 
the semi-infinite ID electrodes, .^/r for the left/right electrodes without the edge atom, 
and J^M/MR for the contact between them. The Hamiltonian for the extended molecule is 
expressed as = €iclci + J2i j ^^jlclc^ -|-h.c.), where is the on-site energy of n orbitals 
for C atoms and of 6s orbitals for Au atoms, is a hopping parameter between ith and 
jth orbitals in equilibrium, and c| (q) is the creation (annihilation) operator of an electron 
on the ith orbital. In the present investigation, the on-site energy eAu for the Au electrode 
is assumed to lie in the center of the HOMO-LUMO gap of the isolated C2o- The hopping 
parameters between Au and C atoms are chosen to be —1.0 eV, whereas those between tt 
orbitals are determined using the TBMD calculations [l^ . 

The electron-vibration interaction is given by 



A ij 



(1) 



where the coupling constant g^j is given as 



4- 



dt 



OL=x,y,z 



2^ 



11^ 



u 



(2) 



where u}^ is an eigenvector of the dynamic matrix and dtij/dsia is the hopping modulation 
with an atomic displacement Sia from equilibrium along the a direction of ith carbon atom. 
The modulations dtij/dsia within the C20 are determined using the TBMD calculations, and 
those between the C20 and Au electrodes are assumed to be 1.0 eV/A. 

According to the NEGF formalism, the electrical current from the left electrode to a 
system is given by 



2f r°° 

I=jJdeTT [S<(6)G>(6) - S>(6)G<(6)] . 



(3) 



For the steady-state transport, the lesser and greater Green's functions obey the steady-state 
Keldysh equation: = {'E^ + 'S^+'Ec-v)G^ , where the self-energy S^^j^^ comes from the 
left (right) contacts and Sllv is related to electron-vibration interaction. The retarded (ad- 
vanced) Green's function is given by the Dyson equation: G^*^"^ = G0''"('*) + G°'''^"^Se-v^G^*^"\ 



where G^''"'-"^ is an unperturbed retarded (advanced) Green's function and Se-v"* is the re- 
tarded (advanced) self-energy due to electron-vibration interaction. In the present calcula- 
tion, the self-energies due to electron-vibration interaction are treated perturbatively using 
the Feynman diagram technique, and are expanded up to the lowest order of self-energy 
diagrams. On the other hand, the self-energies of semi-infinite electrodes can be calculated 
analytically in the case of the semi- infinite ID electrode 18|. The equation © is reduced 
to the well-known Landauer formula for the elastic current I^i by setting 5]^^ = in the 
Keldysh equation. 

The I-V characteristics of the C20 bridges are shown in the inset in Fig. [21 where the 
solid and dashed curves represent the total current / and the elastic current Id, respectively. 
From the calculated I-V characteristics, the total current curve is found to deviate slightly 
upward from the elastic curve as the applied bias-voltage increases. The dissipation power 
into molecular vibrations is also estimated, up to Vuas = 100 mV, as of the order of 1 nW, 
which is less than 3% of the total power generated in the entire bridge. This means that local 
heating due to the molecular vibrations is not a serious bottleneck for transport characteristic 
of the C20 bridge compared with that due to contact resistance. 

The contribution from each vibrational mode to the transport characteristics can be seen 
clearly in the differential conductance, dl/dV^iss, shown in Fig. |21 Solid and dashed curves 
in Fig. 121 represent the total differential conductance dJ/dVbias and its elastic part dld/dVbii^, 
respectively. In contrast to the smooth curve for the elastic differential conductance, large 
discontinuous jumps appear in the total differential conductance curve at particular bias- 
voltages, as indicated by arrows in Fig. |21 The contribution from vibrational modes other 
than those indicated by the arrows is increasingly small, despite the nonzero coupling con- 
stants g^j. A typical example is illustrated by the conductance step at 11.1 mV, originating 
from the shuttle motion of C20, which is two orders of magnitude less than the step at 
29.9 mV. The shuttle motion will be discussed in further detail later. 

The physical origin for the considerable variability of the magnitude of the discontinu- 
ous step in Fig. [21 with the vibrational mode, is now discussed from the viewpoint of the 
electronic structures of the C20 bridge. Fermi's golden rule gives us not only the scat- 
tering rate of electronic states, but also the magnitude of these steps in the conductance 
curve. According to Fermi's golden rule, the magnitude of the steps is proportional to 
I j S'ij-^i ~ ^A)^j(e)P • Here \&i(e) represents a scattering electronic state at the ith 



atom, in the absence of an electron-vibration interaction. The scattering electronic states 
within the narrow bias-window [ep — eV/2,eF + eV/2] can be approximately replaced by 
\l/j(ei7') at the Fermi level, therefore, the magnitude of the conductance steps can be esti- 
mated by Sx = fi'^ \ j 9ijG^j^{ep)\'^ , referred to as the scattering intensity hereafter, where 
^ij^(^i^) — "^^"^^i (fF)^j(eF) is the unperturbed lesser Green's function. Figure |21 shows the 
scattering intensity Sx for the C20 bridge. The intensity exhibits large peaks for particular 
modes indicated by the arrows. These mode energies coincide completely with the positions 
of the discontinuous steps indicated by the arrows in Fig |^ As seen in the expression of 
S\, the magnitude of the conductance steps is determined by the local electronic states near 
the Fermi level, as well as vibrational states of the C20 bridge. G^j'^{eF) describes the cor- 
relation between electrons with the Fermi energy on the ith and jth atoms and the g^j has 
a large value for strong stretching of interatomic bond, therefore, large conductance-steps 
are thought to appear when the bond between atoms with high electron-densities stretch 
strongly. 

Thus, a novel way to control the motion of C20 between two electrodes under current 
flow, utilizing the obtained knowledge from the correlation between scattering electronic 
states and molecular vibrations is proposed. The shuttle motion of C20 is a key motion, 
however, its excitation rate is very low owing to the small Sx at hwx = 11-1 meV in Fig. 01 
As explained previously, the rate is expected to increase significantly when the electronic 
state localized at C atoms adjacent to Au electrodes lie close to the Fermi level, because 
only two springs between C20 and Au electrodes stretch in the shuttle motion. However, 
such a localized state lies 1.5 eV below the Fermi level {ep = eV), as shown in Fig. |31|^a), 
and the scattering intensity of the shuttle motion exhibits a maximum peak at —1.5 eV in 
Fig. life). Therefore, its excitation rate can be enhanced by tuning the gate voltage to shift 
the localized state to the Fermi level. Of course, more precise analyses beyond the rigid-band 
picture assumed here will be necessary to quantify the gate-voltage effect on the Fermi level. 
The control mechanism of the shuttle motion of C20 proposed herein is complementary with 
the Coulomb blockade mechanism for the recent experiment on electronic transports in the 
Cgo weakly suspended between Au electrodes 0]. The mechanism proposed in the current 
investigation is efficient, even for inelastic electronic transport through fuUerenes strongly 
connected to the electrodes. 

We finally emphasize the extention of the proposed idea to other molecular bridges. For 
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the 1,4-benzenedithiol (BDT=C6H4S2) bridge, it has been theoretically predicted that a 
rotational motion of 1,4-BDT molecule is strongly coupled to scattering electronic states at 
low-bias voltages j^^]. By analogy with our results, the rotational motion will be switched 
on when a scattering state with high electron-densities on sulfur atoms in 1,4-BDT adjacent 
to the electrodes is shifted to the Fermi level by tuning the gate voltage. 

In conclusion, the inelastic electronic transport assisted by molecular vibrations in the 
smallest fuUerene C20 bridge was investigated using the NEGF formalism combined with the 
TBMD method. The effect of the molecular vibrations on the transport characteristics of 
C20 appear as discontinuous steps with various heights in the differential conductance curve. 
The stepwise behavior of the conductance, and the variation of magnitude of the step was 
clearly understood by analyzing the electronic structures based on Fermi's golden rule, as 
well as the vibrational structures. Local heating caused by the molecular vibrations of the 
C20 bridge were also found not to be a serious bottleneck for functions of fullerene-based 
nanodevices. Moreover, the idea that the shuttle motion of the fuUerenes can be controlled by 
adequately tuning the gate voltage was proposed. The present study provides a clue toward 
a more complete understanding of inelastic electronic transport and the electromechanics in 
various molecular devices. 

This work was supported in part by "Academic Frontier" Project of MEXT (2005-2010). 
Part of the numerical calculations were performed on the Hitachi SRSOOOs at ISSP, The 
University of Tokyo. 
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FIG. 1: Schematic diagram of the C20 fullerene connected to Au electrodes by two springs. The C20 
with D^d point-group symmetry contains four distinct bond lengths among nonequivalent carbon 
atoms labeled a, 6, and c (or c', c"). 

FIG. 2: The differential conductance in the unit of 2e^//i. Solid and dashed curves represent 
the total differential conductance and its elastic component, respectively. The arrows indicate 
particular discontinuous steps. The inset represents the total and elastic electronic current through 
the C20 bridge. 

FIG. 3: The scattering intensity of electrons due to molecular vibrations of the C20 bridge. The 
arrows indicate the particular peaks that coincide with the discontinuous steps in Fig. |21 

FIG. 4: The electron density of Au-C2o-Au at (a) —1.5 eV and (b) eV (the Fermi level) and 
(c) the scattering intensity for the shuttle motion of C2o- The electron density is indicated by the 
shading on the atom spheres. 
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